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ABSTRACT 

Presented  here  are  results  of  ab  initio  Density  Functional  Theory  (DFT)  structural  relaxations 
performed  on  dehydrogenated  and  monohydrogenated  nanocrystalline  diamond  structures  of 
octahedral  {111}  and  cuboctahedral  morphologies,  up  to  approximately  2nm  in  diameter.  Our 
results  in  this  size  range  show  a transition  of  dehydrogenated  nanodiamond  clusters  into  carbon 
onion-like  structures,  with  preferential  exfoliation  of  the  C( 111)  surfaces,  in  agreement  with 
experimental  observations.  However,  we  have  found  that  this  transition  may  be  prevented  by 
hydrogenation  of  the  surfaces.  Bonding  between  atoms  in  the  surface  layers  of  the  relaxed 
structures,  and  interlayer  bonding  has  been  investigated  using  Wannier  functions. 

INTRODUCTION 

With  the  advent  of  nanotechnology  related  research  in  recent  years,  many  studies  have 
sought  to  understand  the  formation  of  the  carbon  clusters,  and  the  dynamic  phase  changes  and 
stability  relationship  between  graphite  and  nanodiamond  clusters.  It  has  been  observed 
experimentally  [1-3],  that  upon  annealing  nanodiamond  particles  transform  into  onion-like 
carbon,  from  the  surface  inward  (with  the  transformation  temperature  dependant  on  the  size  of 
the  particle).  Such  observations  reveal  preferential  graphitization  and  exfoliation  of  the  diamond 
C(lll)  surfaces  over  other  lower  index  surfaces,  and  the  transformation  of  complete 
dehydrogenated  nanodiamonds  into  onion  carbon.  This  transition  has  also  been  modelled 
theoretically  using  semi-empirical  [2,4-7]  and  ab  initio  methods  [8,9].  Presented  here  are 
structural  relaxations  of  hydrogenated  nanodiamonds  with  octahedral  and  cuboctahedral 
morphologies.  It  is  anticipated  that  saturation  of  the  nanodiamonds  surfaces  with  hydrogen  will 
passivate  the  surfaces  and  stabilize  the  structures  and  preserve  the  diamond  structure. 

AB  INITIO  METHOD 

The  hydrogenated  nanodiamonds  have  been  relaxed  using  the  Vienna  Ab  initio  Simulation 
Package  (VASP)[11,12],  We  used  ultra-soft,  gradient  corrected  Vanderbilt  type  pseudo- 
potentials [13]  as  supplied  by  Kresse  and  Hafner  [14],  and  the  valence  orbitals  are  expanded  in  a 
plane-wave  basis  up  to  a kinetic  energy  cut-off  of  290eV.  The  crystal  relaxations  were  performed 
in  the  framework  of  DFT  within  the  Generalized-Gradient  Approximation  (GGA),  with  the 
exchange-correlation  functional  of  Perdew  and  Wang  (PW91)[15],  A detailed  description  of  this 
technique  may  be  found  in  references  [16,17].  We  have  successfully  applied  this  technique  to 
the  relaxation  of  dehydrogenated  [8-10],  hydrogenated  [18]  and  doped  [10]  nanodiamond 
structures  in  the  past,  with  results  showing  excellent  agreement  with  experiment  and  all-electron 
methods  [19], 
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Figure  1.  The  C35H36,  CimHm  and  CissHino  octahedral  nanodiamond  crystals  after  the  relaxation 
process. 


RELAXATION  RESULTS 

Three  octahedral  nanodiamond  crystals  have  been  generated  by  ‘cleaving’  fragments  from  a 
bulk  diamond  lattice,  and  are  denoted  C35,  C«4  and  Ci65.  These  structures  have  C(1 1 1 )( 1 x 1 ) 
single  dangling  bonds  surfaces,  which  have  then  been  terminated  with  hydrogens  to  produce  the 
C35H36,  C»4H64  and  C165H100  nanodiamonds  (see  Figure  1).  Each  nanodiamond  has  then  been 
relaxed,  and  the  final  structure  compared  to  the  previously  obtained  results  [8.9]  for  the 
dehydrogenated  C35,  Cm  and  C165  nanodiamonds.  Following  relaxation,  the  carbon  framework 
of  the  C35H36  structure  was  found  to  be  almost  identical  to  the  initial  cleaved  fragment,  showing 
no  adoption  of  the  rounded  shape  observed  in  the  relaxed  C35  cluster.  Also,  although  the  Cr4  and 
C165  nanodiamonds  have  previously  been  found  to  transform  into  a two-shell  carbon-onions  [9], 
the  hydrogenated  structures  do  not.  The  CmHm  and  C if,sH  |(Ki  nanodiamonds  retain  the  original 
‘cleaved’  diamond  framework  (with  no  exfoliation  of  the  C(1 1 1)  surfaces),  and  exhibit  little 
change  in  the  length  of  surface  bonds.  The  surface  structure  was  found  to  alter  slightly  during  the 
structural  relaxation,  but  finally  relaxed  to  the  bulk-like  positions.  All  evidence  of  the  phase 
transition  characteristic  of  the  dehydrogenated  [111)  nanodiamonds  has  been  eliminated. 

Two  nanodiamonds  with  cuboctahedral  structure  have  also  been  considered.  These  are  the 
C2oH24and  C142H72  nanodiamonds,  which  are  characterised  by  C(1 1 1 )( 1 xl  ):H  and 
C(100)(2xl):H  surfaces  (see  Figure  2).  Previously  published  results  of  the  C29  nanodiamond 
show  a transformation  of  the  nanocrystal  into  the  C2S  fullerene,  with  an  endohedral  carbon  atom 
(C@C2s)[9].  The  C29H24  nanodiamond  was  however  found  to  be  stable,  and  although  the  bond 
lengths  and  angles  varied  slightly  during  the  relaxation,  the  diamond  lattice  was  preserved  (see 
Figure  2).  Similarly,  the  C|42  nanodiamond  has  been  found  to  exhibit  preferential  graphitization 
of  the  C( 111)  surfaces  [9],  producing  fullerenic  ‘cages’  on  the  surface.  These  fullerenic  cages 
were  observed  briefly  in  the  initial  stages  of  the  relaxation  of  the  C142H72  nanodiamond,  however 
following  a rearrangement  of  the  hydrogens,  the  crystal  relaxed  back  into  the  bulk  diamond-like 
structure  (see  Figure  2). 
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Figure  2,  C29H24and  C142H72  cuboctahedral  nanodiamond  crystals  after  the  relaxation 
process. 


Wannier  functions  and  analysis  of  bonding 

Wannier  functions  provide  a local  and  chemical  view  of  the  bonding,  and  are  related  to  the 
occupied  Kohn-Sham  orbitals  via  a unitary  transformation  [20],  In  periodic  (super-cell)  systems 
this  transformation  is  non-unique  due  to  Bloch  phase  factors  and  Brillouin  Zone  sampling, 
leading  Marzari  and  Vanderbilt  [21]  to  propose  maximally  localised  Wannier  Functions  for 
which  the  total  spatial  spread 
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is  minimized.  In  equation  1,  the  brackets  indicate  expectation  values  with  respect  to  individual 
Wannier  functions  w„(r).  Two  key  pieces  of  information  are  then  available:  the  Wannier 
Function  Center  (WFC)  defined  as  the  expectation  value  of  the  position  operator,  and  the 
localisation  defined  by  the  square  root  of  the  variance  [22],  When  the  Kohn-Sham  orbitals  are 
doubly  occupied  (as  in  this  work),  these  two  quantities  reveal  the  location  and  spatial  extent  of 
each  bond. 

The  Wannier  analysis  was  performed  using  the  Car-Parrinello  Molecular  Dynamics  package 
[23],  taking  as  input  the  relaxed  coordinates  from  VASP.  All  calculations  used  the  BLYP 
exchange-correlation  functional  [24],  Martins-Trouiller  pseudo-potentials  [25],  and  a plane-wave 
basis  set  with  a wavefunction  cut-off  of  476  eV.  Results  for  C29  and  C29H24  are  presented  in 
Figure  3,  showing  the  WFC  as  small  black  circles,  and  the  carbon  and  hydrogen  atoms  as  dark 
and  light  grey  respectively. 

In  the  hydrogenated  structure  the  WFC  are  located  on  the  midline  joining  atom  pairs,  and  so 
indicate  a bond.  For  C-C  bonds  the  WFC  lies  at  the  midpoint,  indicating  a symmetric  charge 
density,  whereas  C-H  bonds  have  more  charge  at  the  hydrogen  end  of  the  bond.  In  either  case 
the  WFC  localisation  is  small  (less  than  1 .7  A),  indicating  that  the  WFC  is  associated  with  only 
two  atoms.  The  Wannier  picture  thus  confirms  o-bonding  with  hydrogenation,  as  expected. 
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Figure  3.  Atoms  and  Wannier  Function  Centers  for  C29  (left)  and  C29H24  (right).  Carbon  and 
hydrogen  are  shown  as  dark  and  light  grey  respectively,  while  small  black  circles  denote  the 
Wannier  Function  Centers. 


A very  different  situation  is  evident  for  the  dehydrogenated  C29  structure,  where  the  absence 
of  hydrogen  leads  to  bent  a-bonds,  distorted  71-bonds  and  an  ionic  central  atom.  Atom  pairs 
connected  via  a single  WFC  indicate  a-bonds,  and  have  localisations  very  similar  to  the 
hydrogenated  case.  However,  none  of  these  WFCs  lie  on  the  midline,  but  instead  are  projected 
radially  outward,  with  the  carbon  atoms  subtending  an  angle  of  approximately  140°.  The 
absence  of  atom  pairs  connected  by  two  WFCs  indicates  that  no  7t-bonding  is  present.  Instead, 
the  inner  WFCs  are  bent  towards  the  central  atom  and  are  significantly  more  delocalised  (~2.6  A 
compared  to  1.7  A for  the  outer  WFCs).  In  a category  of  it’s  own  is  a single  WFC  with  an 
intermediate  localisation  of  1 .9  A,  and  connected  only  to  the  central  atom.  This  indicates  the 
presence  of  ionic  character  in  C29,  in  agreement  with  detailed  ab  initio  calculations  by  Jackson  et 
al  [26]. 


DISCUSSION 

Using  WC  analysis  it  has  been  shown  that  hydrogenation  of  C( 111)  nanodiamond  surfaces 
promotes  a bulk  diamond-like  surface  structure  and  bonding.  As  a preliminary  test  of  the  effects 
of  surface  hydrogenation  on  the  structural  properties  nanodiamond,  the  cohesive  energy  has  been 
considered.  In  a previous  study  [8]  the  cohesive  energy  of  dehydrogenated  nanodiamond  was 
calculated  using  the  total  (spin-polarization  corrected)  energy  per  ion  plotted  against  the 
fractional  ratio  of  dangling  bonds  (Nj/N,),  for  a series  of  dehydrogenated  nanocrystal  structures. 
We  note  that,  for  example,  a dehydrogenated  cuboctahedral  nanodiamond  of  ~8.5nm  diameter 
would  have  Nd/N,  ~ 0.07,  therefore  by  extrapolating  to  the  limit  N(i/N,->0  using  a linear  fit  the 
intercept  gives  an  estimate  of  the  cohesive  energy.  The  cohesive  energy  of  dehydrogenated 
nanodiamond  has  been  calculated  to  be  7.71eV[8]. 
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Figure  4.  Prediction  of  the  cohesive  energy  of  relaxed  hydrogenated  nanodiamond  from  the 
linear  fit  to  the  total  energy  per  atom  versus  the  number  of  hydrogens  per  carbon. 


This  procedure  was  repeated  for  the  relaxed  hydrogenated  nanocrystals,  by  plotting  the  total 
energy  per  carbon  atom  versus  the  number  of  hydrogen  atoms  per  carbon  atom  N|/Nc  as  shown 
in  Figure  4.  By  extrapolating  the  nanocrystal  energies  as  a function  of  the  hydrogen  to  carbon 
ratio  (in  the  limit  as  Nh/Nc->  0)  the  cohesive  energy  for  relaxed,  hydrogenated  nanodiamond  was 
predicted  to  be  7.52eV.  This  value  is  considerably  closer  to  the  calculated  bulk  diamond  value 
of  7.39eV[8],  indicating  that  hydrogenation  of  nanodiamond  also  induces  a more  bulk-like 
cohesive  energy  at  the  nano-scale. 

CONCLUSION 

It  has  been  shown  that  hydrogenation  of  nanodiamonds  C(  111)  surfaces  eliminates  the 
nanodiamond  to  carbon-onion  transition,  and  promotes  the  bulk-diamond  structure  and  bonding. 
This  has  been  illustrated  by  comparing  the  hybridization  of  bonds  present  in  the  relaxed 
hydrogenated  and  dehydrogenated  nanocrystals  via  calculation  of  Wannier  functions.  These 
calculations  confirm  that  while  the  dehydrogenated  structures  contain  distorted  o and  jt-bonds, 
all  of  the  corresponding  hydrogenated  nanocrystals  were  found  to  be  entirely  o-bonded. 
Therefore,  although  it  was  found  previously  that  dehydrogenated  nanodiamonds  with  C(lll) 
faces  were  energetically  unstable,  these  structures  may  be  stabilised  by  saturating  the  surfaces 
following  cleaving.  Hydrogenation  of  the  surfaces  of  the  nanodiamond  structures  investigated 
here  serves  not  only  to  passivate  the  surfaces,  but  also  to  stabilize  the  bulk  diamond-like 
structure,  making  nanodiamond  a suitable  possibility  for  the  nanodevices  of  the  future. 
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